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ABSTRACT 

We report the discovery of a period doubling behaviour in a 2.4 d BL Herculis- 
type variable of the Galactic bulge. Another bulge BL Her-type star (P = 2.25 d) 
is a strong period doubling candidate. Both objects have been identified with the 
OGLE-HI photometry. Possibility of peri od doubling in this t ype o f pulsators has 
been predicted almost twenty years ago bv lBuchler fc Moskalik (1993). Our finding is 
the first observational confirmation of their theoretical results. 

Discovery of the first BL Herculis star showing the period doubling effect motivates 
a new theoretical investigation with the state-of-the-art convective pulsation codes. We 
present the results of initial model survey, specifically aimed at studying the observed 
period-doubled BL Her variable. All of our non-linear models have P = 2 Ad. The 
computations confirm that the period doubling effect is caused by the 3:2 resonance 
between the fundamental mode and the first overtone, as indicated by earlier radiative 
models of Buchler & Moskalik. Comparison of the computed and the observed light 
curves allows to constrain the parameters of the star, in particular its metallicity, which 
appears to be high, Z « 0.01. The recent evolutionary tracks put also constraint on 
the mass of the star, which is close to M = 0.50 M . 

Key words: stars: variables: Cepheids - stars: Population II - stars: oscilla- 
tions - hydrodynamics - stars: individual: BLG184. 7 133264 - stars: individual: 
BLG189. 6 137529 



1 INTRODUCTION 

BL Herculis- type variables are short-period (1 — 4d) sub- 
group of the Type II (or Population II) Cepheids. They 
are predominantly found in the old stellar populations, in 
the globular clusters and in the Galactic halo. They are be- 
lieved to be low-mass stars evolving away from the Zero- 
Age Horizontal Branch (ZAHB) towards the Asymptotic 
Giant Branch (AGB), crossing the insta bility strip (IS) at 
luminosities la rger than RR Lyrae stars (|Strom et al.lll97d : 
iGingoldl [l985l . see also Section [4}. All Type II Cepheids, 
including BL Her variables, obey a well defined period- 
luminosity relation and can be used as distance indicators 
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(e.g. iMaiaess et all l2009al lbl. and references therein). They 
are also key objec ts for studying t he stellar and galactic 
evolution (see e.g.. IWallersteinll2002l 'l. 

Characteristic feature of the BL Her light curves is 
a secondary bump observed in some of the variables. 
The phenomenon is similar to the Hertzsprung bump pro- 
gression observed in the classical, Population I Cepheids 
and is also interpreted as caused by the resonant inter- 
action b etween the fundamental m ode and the second 
overtone l|Hodson. Cox fc Kindll982T ). Some early hydrody- 
namic models of BL Her v ariab les were computed by e 
Hodson. Cox fc Kind (I1982T ) and lCarson fc Stothers! (m: 



Fadevev fc Fokinl ( 1985I - ) computed several Type II Cepheid 



models, including BL Her models. One of their models re- 
vealed oscillations with alternating amplitude and period. 
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A large and systematic s urvey of the radiat i ve BL Her 
models was conduct e d by Buchler fc Moskalikl (1 19921 ) and 
IMoskalik fc Buchlerl dl993ah . In m ost of their model se- 
quences iBuchler fc Moskalikl (|l992h found narrow windows 
of parameters for which the light and radial velocity curves 
exhibited a periodic alternations of deep and shallow min- 
ima - the period doubling effect. The origin of the pe- 
riod doubling in hydrod ynamic models was traced by 
IMoskalik fc Buchlerl (|l990T l to the destabilising role of the 
half- integer resonances, (2n + l)wo = 2ojfc, between the fun- 
damental mode and the higher order overtones. Close to 
the resonance centre period-doubling bifurcation can oc- 
cur: the fundamental mode limit cycle becomes unstable, 
and period-two solution, with alternati ng shape of the light 
and ra dial veloci ty curves can emerge . IBuchler fc Moskalikl 
(|l992h (see also IMoskalik fc Buchlerl Il993bl) showed that 
their BL Her models undergo a period-doubling bifurca- 
tion caused by the 3:2 resonan ce between the fundamen - 
tal mode and the first overtone. IBuchler fc Moskalikl (| 19921 ) 
concluded that period doubling effect, not observed in any 
BL Her star at that time, should be even tually found. More 
recently, iMarconi fc Di Criscienzol (|2007T ) conducted a large 
survey of BL Her-type convective models, however none of 
their models exhibited period doubling. 

Alternating deep and shallow minima are characteris- 
tic feature of RV Tau-type variables, which, together with 
W Vir-type stars, form another subgro up of the Type II 
Cepheids (for recent review see e.g., IWallersteinl 120021: 
ISzabadosl [2010). These longer period variables (periods of 
~ 4 - 20 d for W Vir stars and above ~ 20 d for RV Tau 
stars) are more evolutionary advanced than the BL Her 
variables. The W Vir stars loop into the instability strip 
during the helium-shell flashes, while climbing up along the 
AGB, whereas the RV Tau stars enter the instability strip 
evolving away from the AGB toward a white dwarf domain. 
Period doubling behaviour appears first at P ~ 19 — 20 d; 
this is a borderline between the W Vir and RV Tau classes. 
However, the two classes overlap in period and the distinc- 
tion between them is not clear-cut (e.g., Szabadosl l20ld : 
IWils fc Ote"rol l2008). The RV Tau behaviour was reproduced 
in the radi ative hydrodynamic mode ls of IBuchler fc Kovacsl 

II 19871) and iKovacs fc Buchlerl (|l988l ). As analysed later by 
IMoskalik fc Buchlerl l|l99d ). the period doubling in these 
models is caused by the 5:2 resonance between the funda- 
mental mode and the second overtone. 

Period doubling was also discovered in the Kepler ob- 
servations of RR Lyrae stars exhibiting the Blazh ko phe- 
nomenon l|Kolenberg et al.l l20ld : ISzabo et al l l2010l ). How- 
ever, to the present day period doubling was not detected 
in any BL Her-type star, despite theoretical predictions. 

In this paper we report the discovery of the first BL Her 
variable clearly showing the period doubling phenomenon 
as predicted by the theoretical models. The star was found 
during the preparation of the subsequent part of the OGLE- 

III Catalogue of Variable Stars (OIII-CVS), devoted to the 
Type II Cepheids in the Galactic bulge (Soszyhski et al. 
2011, in prep.). In Section[5]we present and analyse the light 
curve of the star and prove that it is indeed a BL Her-type 
pulsator. We also discuss another candidate in which period 
doubling is strongly suspected, for which however, definite 
claims require more precise photometric observations. 

In Section [3] we present a new set of non-linear hydro- 



dynam ic models. As compared to the Mo skalik fc Buchlerl 
l|l993af ) models, these models are convective and include up 
to date physics. Our efforts focus on modelling the single 
BL Her star with unambiguously detected period doubling. 
In majority of our model sequences we find the period dou- 
bling domains and show that they originate from the reso- 
nant interaction of the fundamental mo de with the first over- 
tone, corroborating the earlier result of lBuchler fc Moskalikl 
(1992). Comparison of the model and observed light curves 
allows to constrain the parameters of the star (mass, lu- 
minosity and metallicity). In Section [4] we confront these 
parameters with the predictions of the stellar evolution the- 
ory. We discuss our results and summarise the conclusions 
in Section [5] 



2 OBSERVATIONS AND DATA ANALYSIS 

In this study we used Cousins /-band photometry of the 
Galactic bulge, obtained during the third phase of the Opti- 
cal Gravitational Lensing Experiment (OGLE-III). The ob- 
servations were collected between 2001 and 2009, with the 
1.3 m Warsaw telescope located at Las Campanas Observa- 
tory, Chile. A limited number of measurements in Johnson's 
V-band was also obtained to determine (V — I) colour of the 
stars. The photometry was reduced with the Difference Im- 
age A nalysis method fDIA. lAlard fc Lupton|[l998l : IWozniakl 
120001 ). Detailed description of the OGLE-III instr umental 
setup and data reduction proced ures can be found in lUdalskil 
|2003l ) and lUdalski et all l|200Sl) . 

During the preparation of the catalogue of Type II 
Cepheids in the Galactic bulge (Soszyhski et al. 2011, in 
prep.), two BL Her stars with unusual properties were iden- 
tified. Both objects display excessive scatter around the light 
curve maxima and minima. This attracted our attention, 
and prompted thorough investigation of these two variables. 



2.1 BLG184. 7 133264: first BL Herculis-type star 
with period doubling 

The OGLE-III /-band photometry of BLG184.7 133264 
(18:02:56.60, -30:42:34.4) consists of 768 measurements, 
spanning 2830 d. We analysed these data using a stan- 
dard consecutive pre-whitening technique. The star displays 
a very pronounced variability with the period of Po = 
2.399267 d (/ = 0.4167939 d _1 ). We fitted this variability 
with the Fourier series of the form: 



I(t) =Io+J2 A k COs(27Tfc/ t + 4> k ) 



(1) 



where the frequency /o was also optimised. In the next step, 
we subtracted the fitted function from the light curve ("pre- 
whitening"), and then searched the residuals for additional 
periodicities. This was done with the standard Fourier power 
spectrum. In Fig. [T] we present the results of this procedure. 
The power spectrum of the pre- whitened data (second panel 
of the figure) is dominated by a peak at f x = 0.20841 d _1 . 
Within accuracy of the data, this frequency is identical 
to 1/2 /o. In other words, f x is a subharmonic of the pri- 
mary pulsation frequency /o. Subsequent step by step pre- 
whitening revealed additional subharmonics, up to 9/2 fo 
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Figure 1. Pre- whitening sequence for BLG184. 7 133264. Up- 
per panel: power spectrum of original data. The spectrum is 
dominated by pulsation frequency fo and its daily aliases. Mid- 
dle panel: power spectrum after removing fo and its harmon- 
ics. Dashed lines indicate locations of the pre-whitened frequen- 
cies. The highest peak corresponds to a subharmonic frequency, 
1/2 /o, its daily aliases are very prominent, too. Two other sub- 
harmonic frequencies (5/2 fo and 7/2 fo) are also visible. Lower 
panel: power spectrum after removing fo, its harmonics and its 
subharmonics. 



(5/2 fo and 7/2 fo are noticeable already in Fig.[T]). No other 
secondary periodicities were found. 

The presence of subharmonic frequencies in the power 
spectrum, i.e. frequencies of the form (n+l/2)/o, is a charac- 
teristic signature of a period doubling. In the time domain it 
means that the light curve of the star is still strictly periodic, 
but it repeats itself not after one, but after two pulsation 
periods. In other words, the light curve displays regular al- 
ternations of higher and lower maxima (and minima). This 
is clearly visible in Fig. O where we plot the light curve 
of BLG184.7 133264 folded with twice the pulsation period 
(2Po)- The difference between the two consecutive maxima 
is small (0.020 mag), nevertheless well visible. Such alternat- 
ing light curves are well known in long er period W Vir an d 
RV Tau variables (P > 19 - 20 d, see lWils fc Oteroll200ct ). 
but they were never before observed in the BL Her stars. 
BLG184.7 133264 is the first BL Her-type variable, in which 
period doubling is unambiguously detected. 

To fully characterise the light curve of 
BLG184.7 133264, we fitted the data with the final Fourier 
series, which included all the detected subharmonics: 

I{t) = Io+J2 Ak cos (27Tfe/ O t + 0fc) + 

k 

+ B n+i/2 cos [27t(n + l/2)/ t + Vn+1/2] • (2) 

n 

This function is plotted in Fig. [51 together with the folded 
photometric data. The amplitude of the main pulsation com- 
ponent is A\ = 0.2742 ± 0.0005 mag. The highest subhar- 
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Figure 2. Light curve of BLG184.7 133264 folded with twice the 
pulsation period, 2Pq- 
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Figure 3. Pre-whitening sequence for BLG189.6 137529. Upper 
panel: power spectrum of original data. The spectrum is dom- 
inated by pulsation frequency fo and its daily aliases. Middle 
panel: power spectrum after removing fo and its harmonics. 
Dashed lines indicate locations of the pre-whitened frequencies. 
The highest peak corresponds to a subharmonic frequency, 1/2 /o, 
its daily aliases are also visible. Lower panel: power spectrum after 
removing fo, its harmonics and its subharmonics. 



monic is over twenty times smaller and has an amplitude of 
B1/2 = 0.0127 ± 0.0005 mag. The power spectrum of resid- 
uals of the fit is displayed in the lower panel of Fig. [T] The 
dispersion of the residuals is a — 0.0086 mag. 



2.2 BLG189.6 137529 : period doubling candidate 

BLG189.6 137529 (18:00:48.31, -29:58:54.4) is another 
BL Her variable in the Galactic bulge, pulsating with the 
period of P = 2.245688 d (/ = 0.4452978 d" 1 ). During al- 
most 8 years of OGLE-III observations the pulsation period 
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Figure 4. Light curve of BLG189.6 137529 folded with twice the 
pulsation period, 2Pq- 

of the star was slowly changing. Therefore, we decided to use 
in the analysis only a subset of data spanning 1540 d, from 
HJD = 2453416 to 2454956 (543 measurements in total). 
The data was analysed in the same way as in the previous 
case. 

The pre-whitening sequence of BLG189.6 137529 is dis- 
played in Fig. [3] Only one secondary peak was detected 
in this star. Its frequency, f x — 0.22257 d _1 , corresponds 
1/2 f Q . Thus, BLG189. 6 137529 shows a signature of a pe- 
riod doubling. This time however, the detected subharmonic 
is extremely weak, with the amplitude of only B1/2 = 
0.0029 ± 0.0004 mag. The light curve of BLG189.6 137529, 
folded with twice the pulsation period is displayed in Fig. [4] 
As we might have expected, the alternations of maxima and 
of minima are barely visible. 

To make them visible better, we display the same data 
again in Fig. [5] but this time we do it differently. We plot 
the folded data twice, with different colours. The two curves 
are shifted by one pulsation period (i.e. by Po). In this way, 
even and odd pulsation cycles are plotted on top of each 
other. The inset provides a closer look at the maxima of the 
folded curves. The red points are systematically higher than 
the black ones, but the difference is very small. The maxima 
of the fit lines differ by only 0.0058 mag, which is less than 
the scatter of the data around the fit (0.0064 mag). 

Both Fig. [3] and Fig. [5] show that the period doubling 
in BLG189.6 137529, is present, but is much less significant 
than in the first star. We consider this BL Her variable to 
be a strong period doubling candidate, which however needs 
further confirmation with additional data. Such data are 
currently being collected (for both stars of this paper) within 
the fourth phase of the OGLE project. 

2.3 Are BLG184.7 133264 and BLG189. 6 137529 
genuine BL Her-type variables ? 

In Fig. |6] we put BLG184.7 133264 and BLG189.6 137529 
on the colour-magnitude diagram, together with almost 
17000 RR Lyr variable s observed toward the Galactic bulge 
jSoszvriski et aljfeoill ). The bulge RR Lyr stars are spread 
over a wide range of (V — I) colours. This is caused by 
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Figure 5. Light curve of BLG189.6 137529 folded with 2P , plot- 
ted twice with different colours. The red and black curves are 
shifted by Pq. The inset shows the difference between even and 
odd pulsation cycles. 



Table 1. Light curves of BLG184.7 133264 and BLG189.6 137529. 
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large field-to-filed variations of interstellar reddening. We 
plot with the filled circles the RRab stars located in the 
same fields as BLG184.7 133264 and BLG189.6 137529. Both 
period-doubled variables discussed in the current paper have 
the same colours as the RRab stars in their respective 
fields (i.e. RRab stars with the same reddening). This im- 
plies that BLG184.7 133264 and BLG189.6 137529 are lo- 
cated inside the instability strip. Their luminosities, which 
are 1.65 ± 0.55 mag brighter than the Horizontal Branch, 
place them in the domain of the BL Her-type pulsators. 
The light curves of both stars are also quite typical 
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Figure 6. Colour-magnitude diagram for BLG184. 7 133264 and 
BLG189.6 137529 (plotted with diamonds). Background grey dots 
represent RR Lyr stars observed toward the Galactic bulge. 
Filled circles mark RRab stars observed in the same fields as 
BLG184.7 133264 and BLG189.6 137529. 
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Figure 7. Light cur ves of selected BL Her stars of the LMC 
(Soszvhski ct al. 2008). Numbers at the right side of the plot show 
periods (in days) of the displayed variables. 



for the BL Her variables. In Fig. we display a selection 
of BL Her stars discov ered in the Large Magellanic Cloud 
|Soszvnski et al.ll2008r ). All these variables were identified 
on the basis of their position in the P — L diagram, which 
established their BL Her nature beyond doubt. Fig. shows 
that there is nothing unusual about the light curves of 
BLG184.7 133264 and BLG189. 6 137529. In fact, they are 
very similar to the light curves of BL Her stars observed in 
the LMC. This statement can be put on a more quantitative 
basis with the help of the Fourier decomposition parameters 
R k i — Ak/Ai and 4> kl — 4> k — k<pi, where the A k and 4>k 
are defined by Eq.[T]or Eq.[2] The values of Rki and <f>ki for 
the light curves of BLG184.7 133264 and BLG189.6 137529 
are given in Tab. [1] In Fig. [8] these values are compared 
with the Fourier parameters of the Type II Cepheids of the 
LMC. Both period-doubled variables are firmly located on 
the Fourier parameter progressions defined by the LMC pul- 
sators. This is particularly striking in the plot of the (f>3i vs. 
period, where Type II Cepheids form a very tight, well de- 
fined sequence. Based on Figs [6] [7] and [8] we conclude that 
BLG184.7 133264 and BLG189.6 137529 are perfectly nor- 
mal BL Her-type variables and, except of the period dou- 
bling, they do not differ in any way from other pulsators of 
this class. 



3 HYDRODYNAMIC MODELS 

The discovery of the first BL Herculis star showing the pe- 
riod doubling motivates a new survey of convective hydro- 
dynamic models including up to date physics. In this Section 
we report the results of such initial model survey in which 
our goals are the following: (i) to check whether the period 
doubling is possible in the convective BL Her models with 
the fundamental mode period equal to 2.4 d; (m) to confirm 
the crucial role of the 3:2 resona nce, 3u;q = 2o>i, evident from 
the earlier radiative models of iBuchler fc Moskalikl (|l992T) : 
(in) to model the light curve of the only firm BL Her star 
showing the period doubling and (iv) to constrain its physi- 
cal parameters (mass, luminosity, metallicity) based on the 
comparison of the computed models with observations. 

In Section 13.11 we briefly present the convective hy- 
drocodes used in this study and provide the details on the 
adopted model parameters. Results of the model survey are 
presented in the following sections, first, the results of exten- 
sive linear model survey (Section 13. 2p . and next, the results 
of non-linear model computations focused on reproducing 
the light curve of BLG184.7 133264 (Section 133)) . 

3.1 Numerical methods 

In all model computation s we use the convect i ve pu lsation 
hydrocodes described by ISmolec fc Moskalikl ((2008) . The 
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Figu re 8. J-band F ourier parameter progressions for the Type II Cepheids. Asterisks represent Type II Cepheids of the LMC 
iSoszvnski et al.ll200gfi . BLG184.7 133264 and BLG189.6 137529 are plotted with diamonds. Their 31 values are shifted up by 2vr. 



codes are Lagrangian. Radiative transfer is treated in the 
diffusion ap proxima t ion. F or the convective energy transfer 
we use the iKuhfurj | 19861 ) time-dependent model adopted 
for the use in the pulsation hydrocodes. Kuhfufi model is 
a simple, one dimensional formulation, with one additional 
equation describing the generation of t he turbulent energy. 
All eq uations and details are given in ISmolec fc Moskalikl 
(2008). The model contains eight order-of-unity scaling pa- 
rameters. These are the mixing-length parameter, a, and 
parameters scaling the turbulent fluxes and the terms that 
drive/damp the turbulent energy, a p (turbulent pressure), 
Q m (eddy- viscous dissipation), a c (convective heat flux), a t 
(kinetic turbulent energy flux), a s (buoyant driving), ay 
(turbulent dissipation) and 7 r (radiative cooling). Theory 
provides no guidance for their values however, some stan- 
dard values are in use. They result from comparison of 
a static, time-independent version of the model with the 
standard mixing-length theory (see IWuchterl fc Feuchtingeil 
1 19981 : ISmolec fc Mosk alik 2008). In practice, values of the 
parameters are adjusted to get the models that satisfy as 
many observational constraints as possible. Parameters used 
in this study are summarised in Table [2] In this initial study 
(see Section [5]) we neglect the effects of turbulent pressure 
and turbulent overshooting. Set PI is our basic choice. 

Static envelope models were constructed using 180 mass 



Table 2. Convective parameters considered in this study. Pa- 
rameters a B , a c , cij, ojp and 7 r are given in the units of standard 
values (a s = a c = 1/2^/273, a d = 8 /3 A /273, a p = 2/3 and 
7r = 2\/3; see lSmolec fc Moskaliisf[2008l . for details). 
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shells, extending down to 2 x 10 K, with the fixed temper- 
ature (T a = 15 000 K) in the anchor zone located 70 zones 
below the surface. Such mesh structure, with large num- 
ber of zones between the hydrogen ionisation region and the 
surface, allows for better resolution of the convective zone 
in the non-linear computations, which becomes very thin 
during the expansion phases, sweeping inward through the 
Lagrangian zones of the model. 

The static model is perturbed ("kicked") with the scaled 
linear velocity eigenvector and time evolution is followed. By 
default, non-linear model integration is conducted for 1000 
pulsation cycles. With typical linear growth rates charac- 
terising BL Her stars, on order of 0.1 d _1 , the limit cycle 
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pulsation is reached only after several tenths of pulsation 
periods. 

Non-linear models of BL Her stars, are numerically 
much more demanding than the models of RR Lyrae stars. 
Having roughly the same mass, BL Her stars are 2-3 times 
more luminous. Extended envelope is strongly non-adiabatic 
and pulsations are rather violent. Consequently artificial vis- 
cosity, for w hich we use the von Neumann and Richtmyer 
prescription (|Richtmverlll957r i. turns on in the models. We 
adopt th e following values for the artificial viscosity parame- 
ters (see lStellingweril[l975l ). C q = 4.0 and a large cut-off pa- 
rameter a cut = 0.1. Large value of the cut-off parameter (0.1 
as compared to 0.0 1 used by default in the r adiative com- 
putations, see e.g. iBuchler fc Moskalikl ll992T ) assures that 
the artificial viscosity plays a role only in the regions of 
the strongest velocity gradients, located in a few outermost 
layers of our models, and consequently does not affect the 
dynamics of pulsatiorQ. 

In order to compare the models with the OGLE /-band 
light curve of BLG184.7 133264, the bolometric light curve 
is transformed to the /-band by computing the bolometric 
correction and (V — I ) colou r at each pulsationphase. To 
this aim static iKurucd {2005) atmosphere modelqj are used. 



3.2 Linear model survey 

In our model survey we consider four values of model masses, 
M = 0.50 M , 0.55 M Q , 0.60 M q and 0.65 M . The adopted 
luminosities are in large range of 90 — 280 L© . In all our 
models the hydrogen abundance is fixed, X = 0.76, and 
three different heavy element abundances ar e considered , 
Z = 0.0001, 0.001 and 0.01. The OP opacities {Seatonll2005l) 
generated for the solar mixture of lAsplund et alJ (|2004l ) 
are used . At the low temperatures th ey are supplemented 
with the I Alexander fc Ferguson ( 1994) opac ity data, which 
is available for the Grevesse fc Noels! ljl993h mixtur^E Pa- 
rameters of the models are not bounded by the evolutionary 
constraints. Later, in Section[4] we compare our models with 
the recent horizontal branch evolutionary tracks. 

Results of the linear model survey are shown in the HR 
diagrams presented in Figs[9land llOI Fig.|9]which shows the 
results for M = 0.55 Mq and Z = 0.001 is more detailed 
than Fig. \W\ in which we show the results for the full grid 
of considered masses and metallicities. In these models con- 
vective parameters of set PI were used. In both figures the 
instability strips are shown as shaded areas, light grey for 
the fundamental mode, and dark grey for the first overtone. 
Typical width of the instability strip is on order of 1 000 K. 
First overtone can be unstable only for the lowest luminosi- 
ties. Thick solid line running across the instability strip is 
the line of constant fundamental mode period equal to 2.4 d 



1 For the models with the highest L/M ratio artificial viscosity 
can contribute to the pulsation damping also in the hydrogen- 
helium ionisation region, but the contribution is small as com- 
pared to the damping caused by the turbulent viscosity. 

2 http://kurucz.harvard.edu/ 

3 The OP and Alexander &; Ferguson opacities are stitched to- 
gether at logT = 3.95, that is in the outer model layers, above 
partial ionization regions driving the pulsation. Consequently the 
mismatch in the abundance mixtures does not affect the model 
behaviour. 
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Figure 9. Results of the linear model survey adopting the con- 
vective parameters of set PI (Table [3J) and for models with 
M = 0.55 M and Z = 0.001. The shaded areas show the in- 
stability strips for the fundamental mode (light grey) and for the 
first overtone (dark grey). The thick solid line is the line of con- 
stant fundamental mode period, Po = 2.4 d, and other lines show 
the loci of various resonances, as indicated in the key. Presented 
results are based on the interpolation in the extensive grid of lin- 
ear models with 10 Lq step in the luminosity and 50 K step in 
the effective temperature. 



(nearly equal to the period of BLG184.7 133264). Thick dot- 
ted line (plotted only in Fig. [9| shows the loci of the 2:1 
resonance 2wo = ^2 which is crucial in shaping the bump- 
progression observed in the light curves of BL Her stars (see 
Section [1]). Other lines show the loci of various half-integer 
resonances, which are most important in the context of pe- 
riod doubling phenomenon. 

In our study we considered all possible half-integer res- 
onances between the fundamental mode and the eight lowest 
order overtones. Loci of all these resonances are plotted in 
Fig. [5] Six resonances cross the Po = 2.4 d line at different 
luminosities and in principle each of them can be responsible 
for the period doubling observed in BLG184.7 133264. The 
most likely candidate is the 3:2 resonance with the first over- 
tone, 3cjo = 2wi. This resonance was proofed to be respon- 
si ble for the period dou b ling in the radiative BL Her models 
of IBuchler fc Moskalikl l| 19921 ). Other resonances however, 
cannot be excluded a priori. Very promising are the high- 
order resonances with the higher-order overtones, with the 
sixth and with the seventh order overtones. These modes 
are special, as they are trapped in the outer envelope, be- 
tween the surface and the hydrogen ionisation front. Ori- 
gin and prope r ties o f thes e modes were studied in de tail by 
IBuchler et all (|l997f ) and IBuchler fc Kollatbl |200ll ). Such 
modes are characterised by a large growth rates as compared 
to their not-trapp ed neighbours. In their stud y of the half- 
integer resonances iMoskalik fc Buchlerl l|l990l ) showed, that 
period doubling behaviour is possible if the resonant mode is 
not damped too strongly (but it does not have to be unsta- 
ble). Consequently, resonances involving the trapped modes 
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Figure 10. The same as Fig. [9] but for a grid of models with different masses (0.50 Mq, 0.55 Mq, 0.60 Mq, and 0.65 Mg , in consecutive 
columns), and metallicities (0.01, 0.001, and 0.0001, in consecutive rows). For clarity, the loci of only three half-integer resonances are 
plotted, as indicated in the key in the bottom, left-most panel. 



are also the likely candidates to cause the period doubling. 
In our models the growth rates of the sixth and seventh 
overtones are close to zero, but the modes remain stable. 
In addition the modes undergo the avoided crossing, which 
may be noticed from the run of res onant lines in Fig. [9] 
(long-dashed lines). We also note that lKollath et all (|201ll ) 
showed that the period doubling discovered rec ently in sev- 
eral R R Lyrae stars showing the Blazhko effect jSzabo et al.l 
|2010T ) is caused by the 9:2 resonance with the ninth order 
overtone, which is a trapped envelope mode in their mod- 
els. Al so our model compu tations for the Blazhko RR Lyrae 
stars (|Smolec et al.l l201ll ) confirm the crucial role of this 
resonance. 

Fig. [10] illustrates how the location of the resonance 
loci vary with mass and metallicity of the models. For clar- 
ity only the three half-integer resonances are plotted in the 
Figure, the 3:2 resonance with the first overtone, 3a>o = 2wi, 
the 7:2 resonance with the fourth overtone, 7ujq = 2o;4, and 
the 9:2 resonance with the seventh overtone (trapped mode), 
9ujo = 2ujr. These three resonances cross the Po = 2.4 d 
line well within the instability strip for almost all consid- 
ered combinations of M and Z. Other resonances often fall 
beyond or close to the edge of the instability strip. It is vis- 



ible that changes in the location of resonant lines are rather 
small, but systematics is clear. Considering e.g. the 3:2 reso- 
nance and its crossing with Po = 2.4 d line, one observe that 
larger the mass (at fixed Z) lower the luminosity and lower 
the effective temperature at the crossing. The same tendency 
is observed with the increasing metallicity (at fixed mass). 
The largest possible luminosity at the crossing is L « 185 L0 
at T cff « 6316 K (M = 0.50 M , Z = 0.0001), while the 
lowest possible luminosity is L m 147 L© at T e g ~ 5804 K 
(M = 0.65 M , Z = 0.01). We note that model periods de- 
pend only weakly on the convective parameters of the mod- 
els. Sensitivity is larger for the period ratios and thus, also 
for the location of the resonant lines. The quantitative pic- 
ture arising from Figs [9] and [lJJ however, is the same for 
other sets of convective parameters, in particular for sets 
P2 and P3 considered in this study. 

Linear computations show, that at least several half- 
integer resonances can play a role in destabilising the 
fundamental mode and giving rise to the period doubling 
behaviour in a models with fundamental mode period 
of 2.4 d. However proximity to the resonance centre and 
weak damping of the resonant overtone are not sufficient 
conditions. The non-linear effects which determine the 
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mode interaction are crucial here and only full amplitude 
non-linear computations can reveal whether and which 
resonance can cause the period doubling. 



3.3 Non-linear model survey 

In this Section we present the results of non-linear model 
survey aimed to model the light curve of BLG184.7 133264, 
in particular, to reproduce the period doubling behaviour 
and to investigate its causes. All model sequences run across 
the fundamental mode instability strip, along the lines of 
constant period equal to 2.4 d (Figs 191 and [ID) . Thirty six 
model sequences were computed, for all considered values 
of model masses (0.50 M Q , 0.55 M , 0.60 M and 0.65 M ), 
model metallicities (0.01, 0.001 and 0.0001) and for the three 
sets of convective parameters (PI, P2 and P3, Table [2). In- 
dividual models in a sequence are specified by the static 
luminosity, which is changed in the 2Lq steps. Effective 
temperature is adjusted to match the Po = 2.4 d (for con- 
sidered mass and metallicity) . We note that full amplitude 
non-linear periods may slightly differ from the linear value 
and in our models are usually longer by up to 3 per cent. 

Bolometric light curves are transformed to the /-band 
(see Section l3.ip and next, are fitted with the Fourier se- 
ries (equation [2) . We also compute the Fourier decomposi- 
tion parameters of low order, the amplitude ratios, -Rfei, and 
phase differences, <pki (see Section f2 . 3 [1 . 

Exemplary /-band light curves from model sequence 
with M = 0.50 M Q and Z = 0.01 (set PI) are shown in 
Fig. 1111 Pulsation period is roughly the same for all these 
models, Po ~ 2.4 d. Systematic change of the light-curve 
shapes is clearly visible. Model with the largest luminosity 
(190 Lq) has a regular triangular shape almost void of sec- 
ondary features. As luminosity decreases the secondary fea- 
tures such as bumps, become more and more pronounced. Of 
great interest is the bump close to the minimum brightness, 
which is also visible in the light curve of BLG184.7 133264 
(Fig. [2J), as well as in some other BL Her- type variables 
(e.g. Fig. [T) . Its amplitude increases with the decreasing lu- 
minosity. Another bump is visible close to the maximum 
brightness and is most pronounced for the lowest luminosity 
models. One may also note a series of wiggles of very low 
amplitude on the descending branch. It is a numerical ef- 
fect caused by the lack of resolution in the convective zone 
which becomes very thin during the expansion phases (close 
to maximum radius, see also Section l"3.1[) . Period doubling is 
clearly visible in several models with luminosities in a range 
140 - 170 L Q . 

For other model sequences we observe qualitatively the 
same picture. In nearly all our model sequences period dou- 
bling behaviour is evident for several models. In Section l"3.3,ll 
we focus on the origin of this phenomenon. Later in Sec- 
tion [310] we focus on the light-curve shapes and search for 
the best-fitting model for BLG184.7 133264. 

3.3.1 Origin of the period doubling 

Period doubling is caused by the destabilisation of the fun- 
damental mode by the resonant overtone. The stability of 
the limit cycle pulsation can be s tudied using the relax- 
ation technique (|StellingwerJll974l ), which shows whether 
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Figure 11. A sample of /-band light curves from the model se- 
quence with M = O.5OM and Z = 0.01 (set PI). For clarity, 
the consecutive light curves are shifted arbitrarily by 0.3 mag. 
Each light curve is marked with luminosity value on the left- 
hand side. Light curves showing the period doubling phenomenon 
(140 — 170 Lq) are plotted with thick lines. 



the non-linear limit cycle is stable or not with respect to 
the possible perturba tion in other (resonant) m ode. Such 
analysis was done by iBuchler fc Moskalild ij 19921 ) for their 
radiative BL Her models. The period doubling domain they 
found coincide with the 3:2 resonance centre. Using the re- 
laxation technique they showed that for these models the 
fundamental mode limit cycle is unstable with respect to 
the first overtone, pointing to the crucial role played by the 
resonance. 

Unfortunately, the relaxation scheme is not imple- 
mented in our convective hydrocodes. Consequently, we are 
left to judge on the role played by particular half- integer res- 
onance in the computed period doubling behaviour, based 
on the distance of the model to the resonance centre (mis- 
match parameter, A), which is computed using the linear 
period ratios. One have to remember that when non-linear 
results are correlated with the linear properties of the mod- 
els, the models have to be computed in a consistent way, i.e. 
using th e same numerical m esh and the same input physics 
(see e.g. IBuchler et aLlI 19901 1 . This is the case of our models. 

The existence and strength of the period doubling be- 
haviour in our models is characterised by the amplitude of 
the sub-harmonic frequency, B 1 / 2 (eq. [2). In Fig. [12] /2 is 
plotted versus the mismatch parameter for the three half- 
integer resonances, 3ljo — 2wi (A3 : 2 = cji/cjo — 1.5) in the 
top panel, 7luq — 2oj4 (A-7-,2 = oja/ojq — 3.5) in the middle 
panel, and 9u;o = 2lu7 (Ag : 2 = ujt/ojo — 4.5) in the bottom 
panel. Each line in this figure corresponds to one model se- 
quence with given mass and metallicity. All model sequences 
displayed in Fig. 1121 adopt convective parameters of set PI. 
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Figure 12. Amplitude of the sub-harmonic frequency compo- 
nent, plotted versus the mismatch parameters for: (top 
panel) 3:2 resonance with the first overtone, A3 : 2 = ujx/uiq — 1.5, 
(middle panel) 7:2 resonance with the fourth overtone, Ay : 2 = 
014/0)0 — 3.5 and, (bottom panel) 9:2 resonance with the seventh 
overtone, Ag ; 2 = 0)7/0)0 — 4.5. Models with different masses are 
plotted with different line styles as indicated in the key. Metal- 
licity is indicated with the line thickness; the thick lines are for 
the largest metallicity, Z = 0.01 and the thin lines are for the 
lowest metallicity, Z = 0.0001. For clarity, model sequences with 
Z = 0.001 are not shown in the Figure. In all models convective 
parameters of set PI (Table [2} were adopted. 



For all model sequences we see a single period doubling do- 
main with typical amplitudes of few hundreds of magnitude. 

In the top panel of Fig. [12] we clearly see that the pe- 
riod doubling domains always coincide with the 3ojo = 2wi 
resonance centre. In most cases, the resonance centre is 
located within the period doubling domain. For all mod- 
els showing the period doubling, the mismatch parameter 
is always smaller than 0.05 (and is usually much smaller). 
We note that period doubling domain can be significantly 
shifted with respect to the resonance centre, as it is e.g. for 
M = 0.60 M and M = 0.65 M (Z = 0.01) domains which 
lie outside the resonance centre. It does not disproof the role 
played by the resonance however, but can be caused by the 



non-linearity and non-adiabaticity of the pulsations (see e.g. 
iBuchler fc Moskalik|[l992l : ISmoledl2009al ). 

The period doubling domains are also close to the 
7luq = 2u)4 resonance centre (middle panel of Fig. I12[) . In 
most cases the resonance centre falls somewhere within the 
period doubling domain. But contrary to the 3wo = 2lji 
resonance, the range of mismatch parameters is very large, 
—0.5 < A7;2 < 0.2. It is highly unlikely that the period dou- 
bling is caused by the 7ojo = 2oj4 resonance in the model 
located as far from the res onance centre as e.g. | A7 : 2 1 > 0.2 
|Moskalik fc Buchlerlfl990l ). 

For the 9ljo = 2oji resonance (bottom panel of Fig. I12[) 
conclusions are the same. Most of the period doubling mod- 
els are located far from the resonance centre. Even the entire 
period doubling domain can be located as far as Ag : 2 > 0.15 
(for model sequence with M — 0.50 Mq and Z — 0.01 - thick 
solid line in the bottom panel of Fig. I12[) . Clearly, the 9:2 
resonance cannot play any role in these models. 

Similar analysis was done for all other possible half- 
integer resonances which centres fall within the computed 
model sequences (e.g. for the other three resonances which 
loci are plotted in Fig. [9]), as well as for all model sequences 
adopting the convective parameters of sets P2 and P3. Only 
for the 3:2 resonance the mismatch parameter is very small 
(1^3:2 1 < 0.05) for all the models that show the period dou- 
bling. It clearly indicates that the 3:2 resonance plays a 
key role in causing the phenomenon. The question whether 
other resonances can play an additional role for some par- 
ticular models, for which respective mismatch parameter is 
also close to 0, cannot be answered beyond doubt. It seems 
however, that this is not the case. The period doubling do- 
mains as depicted in Fig. [12] through the amplitude of the 
sub-harmonic frequency, -B1/2, seem to be void of features 
correlated with the location of the centres of the other half- 
integer resonances (see also Fig. I13[) . The crucial role of the 
3:2 resonance in causing the period doubling behaviour in 
BL Her type models is al so supported by the earlier analy- 
sis of radiative models by Buchlcr & Moskalik (1992). Using 
the relaxation technique lBuchler fc Moskalik! (|l992T ) proofed 
the key role of the 3:2 resonance. 

In our opinion, the presented results leave no doubt that 
the 3:2 resonance with the first overtone is responsible for 
the period doubling in the hydrodynamic models and first 
of all in BLG184.7 133264. 



3.3.2 The best model for BLGI84. 7 133264 

In this Section we compare the computed model light curves 
with the light curve of BLG184.7 133264 (Fig. O . Our com- 
parison is based on the low-order Fourier decomposition pa- 
rameters, the amplitudes, A\ and -B1/2, the amplitude ratio, 
R21, and the phase difference, 4>2\. Higher order parameters 
are modelled less accurately with the hydrocodes and are 
not used below. We note that physical parameters of all 
our models were adjusted to match the linear fundamental 
mode period of 2.4 d with typical accuracy larger than 0.01 d. 
Non-linear periods may differ only very slightly. These pe- 
riods are close enough to the period of BLG184.7 133264 
and thus allow the direct comparison of the Fourier pa- 
rameters which is done in Fig. Q3] {A\ versus -B1/2) and in 
Fig. [14] (021 versus i?2i). In these figures we show the re- 
sults for all the computed models, for the full grid of model 
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Figure 13. Amplitude of the sub-harmonic frequency component, B1/2, plotted versus the amplitude, A±, for the full grid of the 
computed non-linear models. Different sets of convective parameters (Table [2} were used in the computations, as indicated in the key 
in the top, left-most panel. Cross shows the amplitudes of BLG184. 7 133264 (see Section [2). The best matching models of set PI, A 
(M = 0.50 M Q , Z = 0.01), B (M = 0.55 Mq, Z = 0.01) and C (M = 0.60 M©, Z = 0.01), are encircled and labelled. 



masses (0.50 M , 0.55 M , 0.60 M Q and 0.65 M Q , consec- 
utive columns of Figs [13] and I14[) and model metallicities 
(0.01, 0.001, and 0.0001, consecutive rows of Figs \TS\ and 
I14[) . Model sequences adopting different sets of convective 
parameters are plotted with different symbols, circles (set 
PI), squares (set P2) and triangles (set P3). Filled sym- 
bols correspond to the models showing the period-doubling 
behaviour, while the models represented by open symbols 
show no sign of alternations. Only the former models are 
shown in Fig. 1131 By default, the consecutive models in each 
model sequence differ by 2 Lq . Sometimes the difference is 
larger, as some models could not be converged to the limit 
cycle due to numerical problems (this is true particularly 
for set P3). Cross in the figures shows the parameters of 
BLG184.7 133264 (Tab.p. 

To find the best matching model for BLG184.7 133264 
we focus on the model amplitudes (Fig. I13|) . the pulsation 
amplitude, Ai, and the amplitude of period doubling, -B1/2. 
One could choose the Fourier phases instead (<f>2i and (fizi), 
which little depend on the pulsation amplitudes and are 
also n ot so sensitive to the values of convective param- 
eters <|Smoledl2009bh . We note however, that one of the 
known problems in modelling the light curves with pulsa- 
tion hydrocodes is that computed Fourier phases are sys- 
tematically different from the observed phases. The prob- 
lem was studied in more deta il in the RR Lyrae stars 
(e.g. | Dorfi fc Feuchtingerl 1 19991 ; iKovacs fc KanbuJ 1 19981 ; 
iFeuchtingerl 1 19991 ). It was also encountered in the ra- 



diative and convective models o f classical Cep heids (e.g. 
iMoskalik, Buchler fc Maromlll992l;ISmoledl2009bl) and in ra - 
diative models of BL Her stars (|Moskalik fc Buchleilll993an . 
Origin of the discrepancy is not clear, although the very 
simplified treatment of the radiation transfer in the subpho- 
tospher ic layers of all current hyd rocodes is the most likely 
culprit (|Dorfi fc Feuchtingenll999T ). 

Consequently we choose to focus on the amplitudes first, 
even though they depend on the convective parameters of 
the models stronger, and in the next step, we compare the 
Fourier phases and amplitude ratios. 

In Fig. [T3] we observe that in almost all model se- 
quences the period doubling domain is present. Only for 
three model sequences, all with the largest considered metal- 
licity (Z = 0.01, and M = 0.55 - 0.65 M©), and adopt- 
ing convective parameters of set P2, period doubling is not 
present. The other striking feature visible in Fig. [T3] are 
large amplitudes of the period doubling as compared to 
BLG184.7 133264 and, at the same time, significantly lower 
pulsation amplitudes, A%. The systematics is clear. Lower 
the mass and higher the metallicity of the models, higher 
the pulsation amplitude and lower the amplitude of period 
doubling. Only for the largest metallicity (Z — 0.01) we can 
match the pulsation amplitude of BLG184.7 133264 (for sets 
PI and P3). Also for the most metal rich model sequences 
the amplitude of period doubling is comparable to what we 
observe in BLG184.7 133264. For lower metallicity models 
(Z = 0.001 and Z = 0.0001) the period doubling amplitude 
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Figure 14. The phase difference, <f>2\ , plotted versus the amplitude ratio, R21 , for the full grid of the computed non-linear models. 
Different sets of convective parameters (TablefS} were used in the computations, as indicated in the key in the bottom, left-most panel. 
Filled symbols correspond to the models showing the period doubling. Crosses represent the Fourier parameters of BLG184.7 133264 (see 
Tablc[T]l. The best matching models of set PI, A, B and C, are encircled and labelled. 



is much higher through the whole period doubling domain 
and matches the amplitude of BLG184.7 133264 only at the 
domain edges. The large amplitude of period doubling in our 
models and rather wide range of luminosities over which 
period doubling domains extend (typically 30 — 40 Lq) is 
somewhat surprising. This suggests that the period doubling 
should be a rather common feature of BL Her stars with pe- 
riods around 2.4 d, which is not the case. We will address 
this problem later in Section [5] 

The above analysis of model amplitudes points that 
BLG184.7 133264 should be of high metallicity. Then, we 
can simultaneously match its pulsation amplitude and the 
amplitude of period doubling. This is a likely possibility, as 
many BL Her stars have high metallicities (see discussion 
in Section [5}. However, one have to remember that model 
amplitudes depend on the convective parameters adopted 
in the computations, in particular on the mixing length and 
on the eddy-viscosity parameter. Comparison of the model 
results for sets PI and P2 reveals the effect of changing the 
mixing length parameter (while keeping all other convective 
parameters fixed). Mixing length is equal to 1.5 in our basic 
set PI and is larger in set P2 (a — 1.8). The effect is clear, 
by increasing the mixing length the pulsation amplitude de- 
creases as well as the amplitude of period doubling (which 
vanish for Z = 0.01 and three largest considered masses, 
Fig. I13p . It seems that one can match the pulsation ampli- 
tude of BLG184.7 133264 also for lower metallicities by the 
decrease of the mixing length below the value adopted in set 



PI (a — 1.5). However at the same time the amplitude of 
period doubling, which is already very large, would increase 
even more. Also values of the mixing length below 1.5 seem 
too small, typical values used in model computations being 
around 1.5—1 .8. Recent numerical simula tions point to even 
larger values l|Trampedach fc SteinllioTlh . 

Mixing length parameter is usually fixed in model com- 
putations and it is the eddy viscosity parameter, which is 
commonly used to adjust the pulsation amplitudes in the 
convective models. It regulates the amount of viscous dis- 
sipation in the model and is a physically motivated coun- 
terpart of the artificial viscosity in the radiative models. 
Pulsation amplitudes of RR Lyrae and Cepheid model s 
clearly depend on that parameter (see e.g. ISmoledl2009bh . 
Lower the eddy-viscous dissipation larger the pulsation am- 
plitude. Also the existence and the amplitude of the period 
doubling computed in the convective RR Lyrae models by 
iKollath et all (|201lh depend on the eddy viscosity parame- 
ter. Again, lower the eddy-viscous dissipation larger the am- 
plitude of the period doubling in case of the RR Lyrae mod- 
els. To match the pulsation amplitude of BLG184.7 133264 
for metallicities lower than Z — 0.01 one has to decrease 
the eddy viscosity, which will also increase the amplitude of 
period doubling. 

To test the effect of decreasing the eddy-viscosity pa- 
rameter we have computed two additional model sequences, 
both with M = 0.55 M Q and Z = 0.0001, and with eddy- 
viscosity parameter equal to a m = 0.20 and a m = 0.15. The 
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Figure 15. Amplitude of the sub-harmonic frequency compo- 
nent, -B1/2, plotted versus the amplitude, A\, for the three model 
sequences, all with M = 0.55 Mq and Z = 0.0001, adopting dif- 
ferent values of the eddy viscosity parameter: a m = 0.25 (solid 
line), a m = 0.20 (dashed line) and o m = 0.15 (dotted line). All 
other parameters of the turbulent convection model are the same 
as in set PI. 

results are shown in Fig.[l5]in which B 1 / 2 is plotted against 
Ai for models of set PI (solid line, a m = 0.25) and for two 
additional model sequences, with eddy-viscosity parameter 
decreased to 0.20 (dashed lines) and 0.15 (dotted line). We 
note that further decrease of eddy viscosity is not easy. Al- 
ready with a m = 0.20 pulsation become very violent and 
we start to face the convergence difficulties. It is clear that 
amplitude of alternations strongly increase with the decreas- 
ing eddy viscosity parameter. However pulsation amplitude 
increased only very weakly - much too small to match the 
amplitude of BLG184.7 133264. 

Changes of the other convective parameters also affect 
the pulsation amplitudes which is usually compensated by 
the change of the eddy viscosity parameter. This is what we 
have done in set P3, in which we consider the radiative cool- 
ing of the convective elements. Inclusion of radiative cooling 
makes the convection less efficient. Consequently the pulsa- 
tion amplitudes increase which has to be compensated with 
the larger value of the eddy viscosity parameter. For the 
model sequences adopting convective parameters of set P3 
we get the period doubling amplitudes higher than for set 
P2 and similar to that of set PI. The pulsation amplitudes, 
A\, are similar as in set P2. Set P3 is interesting because it 
leads to significantly different light curve shapes than sets 
PI and P2, which we discuss in the next paragraphs. 

Based on the model amplitudes we select three models, 
which match the amplitudes of BLG184.7 133264, for the 
further, more detailed comparison. The models, A, B and 
C, are encircled and labelled in Fig. 1131 All models adopt 
convective parameters of set PI, have Z — 0.01 and dif- 
fer in mas (0.50 Mq for model A, 0.55 Mq for model B and 
0.60 Mq for model C). Physical parameters of these mod- 
els are collected in the upper section of Table Before 
discussing these models in more detail, we turn into the 
analysis of the low order Fourier decomposition parameters, 
the R21 amplitude ratio and the $21 phase difference. In 
Fig. [14] the two quantities are plotted for the full grid of the 
computed models. The range of variation of R21 and (f>2i is 
very large for the considered models all of which have the 
same pulsation period (« 2.4 d). Relation between the two 
Fourier parameters is very similar for sets PI and P2, which 



differ in the value of the mixing length only, and is signifi- 
cantly different for set P3 in which the effects of radiative 
losses are turned on. The period doubling domains (filled 
symbols in Fig. I14[) extend over a large range of Fourier 
phase, (j>2i- There are period doubling models with Fourier 
phase below 5.0, as well as with Fourier phase well above 
6.5. Consequently in all our model sequences we can point 
the model which shows the period doubling and for which 
4>2i matches the value of BLG184.7 133264. For the ampli- 
tude ratio we observe that period doubling domains extend 
around the minimum of R21 in the given model sequence. 
The simultaneous match of the BLG184.7 133264 values of 
R21 and 4>2i (cross in Fig. I14[> is possible for several model 
sequences, which are located along the diagonal in the mass- 
metallicity plane displayed in Fig. 1141 The best match for 
the lowest metallicity sequences (Z = 0.0001) is possible for 
the highest masses (0.60 Mq and 0.65 Mq). For intermediate 
metallicity (Z = 0.001) we get the best match for the inter- 
mediate masses (0.55 Mq and 0.60 Mq), and for the highest 
metallicity (Z = 0.01) we get the best match for the lowest 
masses (0.50 M ). 

Although for several models we can get an almost exact 
match between the observed and model values of both R21 
and 02i parameters, only for the highest metallicity mod- 
els we can match the amplitudes of BLG184. 7 133264. For 
lower metallicities (Z = 0.001 and Z = 0.0001) the pulsa- 
tion amplitudes are lower than observed and the amplitudes 
of alternations are much higher (see Fig. 1 13 p . 

The three models of Z — 0.01, A, B and C, which match 
the amplitudes of BLG184.7 133264 best, are also labelled in 
Fig. 1141 These models are relatively distant from the cross 
showing the location of BLG184.7 133264 and for the lower 
luminosity models from the same model sequence the -R21 
and (j>2i can be matched much better. Still the agreement for 
models A, B and C is reasonable, as quantitative comparison 
presented in the lower section of Table [3] shows. For model 
C the -R21 amplitude ratio is lower than observed by only 
5 per cent and the 021 phase difference is lower by 0.41 rad. 
The agreement is similar for model B and slightly worse for 
model A for which the amplitude ratio is significantly differ- 
ent than observed. In Table [3] we also provide the comparison 
for higher order Fourier parameters, R31 and </>3i. The dif- 
ferences are very large. Model values of R31 are much larger 
than observed, even for the best model C, by as much as 
85 per cent. For 1^31 the lack of agreement is also remark- 
able, the model phases being much lower than observed. 
Detailed analysis of all our models points that although we 
can easily find models for which two Fourier parameters (e.g. 
both phases) simultaneously match the observed values, the 
simultaneous match of all low order Fourier parameters (am- 
plitudes, amplitude ratios and phases) is not possible. 

In Fig.[l6]we directly compare the /-band light curves of 
BLG184.7 133264 (Fig.^ and the three selected models A, B 
and C. The model light curves were shifted vertically in order 
to match the maximum brightness of BLG184.7 133264, and 
shifted in phase, so the ascending branches of the light curves 
overlap. The overall light curve shape of BLG184.7 133264 
is nicely reproduced by all three models, but several dis- 
crepancies are apparent. In the observed light curve, after 
the pronounced shoulder on the ascending branch the light 
raise becomes less steep. The observed light curve is rounded 
at the maximum light. The shoulder is also visible in the 



14 R. Smolec et al. 



Table 3. Parameters of the best matching models for 
BLG184.7 133264. In the top section of the Table physical pa- 
rameters are provided and in the bottom section of the Table the 
Fourier decomposition parameters are given. The differences with 
respect to BLG184.7 133264 (model— star) are given in the paren- 
theses (relative differences for amplitudes and amplitude ratios). 



model: 


A B 


C 


physical parameters 


M/Mq 


0.50 0.55 


0.60 


L/L e 


168.0 176.0 


184.0 


T eS [K] 


6241.6 6201.4 


6167.2 


Fourier parameters 


A 1 


0.2784 (+1.5%) 0.2740 (-0.1%) 


0.2667 (-2.7%) 


Bl/2 


0.0135 (+6.5%) 0.0131 (+2.8%) 


0.0103 (-19.1%) 


flai 


0.304 (+40.7%) 0.245 (+13.4%) 


0.205 (-5.1%) 


4>21 


5.010 (-0.451) 5.022 (-0.439) 


5.060 (-0.401) 


R31 


0.236 (+211%) 0.185 (+143%) 


0.140 (+84%) 


031 


3.376 (-0.618) 3.054 (-0.940) 


2.709 (-1.285) 



models however, the ascending branch is steep all the time 
till the maximum brightness which consequently occurs at 
earlier phases in the models. The light fall also begins ear- 
lier and as a result the overall width of the light curve is 
slightly smaller in the models. At brightness minimum the 
pronounced bump is clearly visible in both the observations 
and in the models. Its location is well reproduced by the 
models, but its amplitude in the models is much higher. 
The changing character of the brightness minima around the 
bump during the alternations (each second minimum after 
the bump is deeper) is qualitatively reproduced by the mod- 
els. Clearly we get the best agreement for model C, which is 
also confirmed by more quantitative comparison (Table [3]), 
but the agreement for the other two models is reasonable, 
as well. 

The discrepancies between the model and the observed 
light curves are not very serious. They concern only the 
more subtle features of the light curve, like the detailed 
shape of the brightness maximum or minimum. Such dis- 
crepancies are reflected predominantly in the higher order 
Fourier parameters. Accurate modeling of these higher or- 
der Fourier terms is a very challenging problem, though. Of 
particular importance is the treatment of radiative transfer 
in the outer zones of the model and the transformation of 
the bolometric light curve to the /-band light curve. These 
are both treated in a very approximate way in our code (see 
Section |3.1[) . Therefore, we do not expect that the models 
will satisfactorily reproduce all details of the observed light 
curves. This is a well known deficiency, common to all hy- 
drodynamic models of la rge amplitude radial pulsators (e.g., 
iKovacs fc Kanburlll99Sl ). 

Comparison of the observed and model light curves 
for BLG184.7 133264 point that the star should be of rel- 
atively high metallicity [Z ~ 0.01) and its mass should be 
lower than 0.60 Mq. Its luminosity should be large around 
170 — 180 Lq. In the next Section we confront these parame- 
ters with the results of recent horizontal branch evolutionary 
models. 
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Figure 16. Comparison of the observed /-band light curve of 
BLG184.7 133264 (solid line) with the best matching models of 
0.50 M Q (top panel), 0.55 Mq (middle panel) and 0.60 M Q (bot- 
tom panel). The model light curves were shifted vertically to 
match the maximum brightness of BLG184. 7 133264 and shifted 
in phase so the ascending branch of the models and observed light 
curve overlap. 



4 COMPARISON WITH THE 
EVOLUTIONARY TRACKS 

The evolutionary status of BL H er stars was di scussed in a 
series of papers by Gingold (see lGingoldlll985l , for a sum- 
mary). Since then, several groups published the horizon- 
tal branch (HB) evolutionary tracks for extensive grid of 
models covering a large range of masses a nd chemical com - 
positions. These i nclude e.g. Y 2 tracks by lYi et al.l (|l997l ). 
Pa dova tracks by iGirardi et al.l (|2000h . the BaSTI tracks 
by IPietrinferni et all l|2006l ) and the Dartmouth tracks by 
iDotter et alj i|200Sl) . Below we briefly review the Gingold's 



results and compare them with the outcome of more recent 
evolutionary computations. To this aim we use the BaSTI 
and Dartmouth tracks. Based on the location of the BL Her 
instability strip in the HR diagram we conclude on the pos- 
sible masses of BL Her stars as predicted by the evolution 
theory. Finally, a comparison with the parameters of our 
best fitting models for BLG184.7 133264 (Section |3~X2"|| is 
done. 

Modelling of the post-main sequence evo lution of the 
low-m ass stars is particularly difficult (see e.g. IChiosi et al.l 
1992). During the evolution along the red giant branch 
(RGB) these stars develop a degenerate helium core. Then 
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mass of the core reaches a critical value the runaway helium 
burning (so called helium flash) is initiated. This violent and 
rapid phase of the evolution, accompanied by a substantial 
mass-loss, cannot be followed with the current evolutionary 
codes. Instead the evolution is restarted at the Zero-Age 
Horizontal Branch (ZAHB), at much larger effective temper- 
atures. Progenitors of the BL Her stars start their evolution 
at the blue side of the ZAHB, far beyond the blue edge of the 
instability strip. During the core helium burning, which is a 
relatively long phase of the evolution, they remain close to 
the ZAHB. As the helium becomes gradually depleted in the 
core, the evolution becomes faster and star moves towards 
the red, crossing the instability strip and thus, becoming a 
BL Her variable. According to Gingold evolutionary compu- 
tations two scenarios are possible, w hich depend m ostly on 
the mass of the envelope (see fig. 1 in lGingoldlll985t ). In the 
first scenario, the star crosses the instability strip once and 
evolves away up along the asymptotic giant branch (AGB). 
This is the scenario which one can also see in Figs ll7l and ll8l 
In the second scenario, so called 'bluenose' is possible. Due 
to the structure changes between the helium and hydrogen 
burning shells, after approaching the AGB the evolution is 
directed to the blue and then, back towards the AGB. Two 
additional crossings of the instability strip are then possible 
before the star finally climbs-up along the AGB. In both 
scenarios the star can enter the instability strip again, dur- 
ing the post-AGB evolution (while executing the thermal 
pulses), becoming a W Vir or RV Tau-type variable. 

Both types of tracks can be found in the more re- 
cent evolutionary computations, however, the 'bluenose' be- 
haviour is very rare. It occurs in a narrow range of masses, 
and only for larger than solar metallicities (BaSTI, Dart- 
mouth) and in addition in models assuming a large enhance- 
ment of the a-elements ([a/Fe] > +0.6, Dartmouth). We do 
not discuss this scenario in the following. 

The recent evolutionary tracks are displayed in Figs 1171 
(BaSTI) and 1181 (Dartmouth ). In these figures the edges of 
the instability strip are marked with solid, vertically run- 
ning lines. The grey-shaded area in between corresponds to 
the BL Her pulsation domain with fundamental mode peri- 
ods between 1 d (lower boundary of the shaded area) and 4 d 
(upper boundary of the shaded area)Q. Below that region the 
RR Lyrae instability strip extends. The following discussion 
is based on the BaSTI tracks. In the three panels of Fig. [T7] 
we have plotted the evolutionary tacks for three different 
metallicities, Z = 0.0001 in the left-hand panel, Z = 0.001 
in the middle panel and Z — 0.01 in the right-hand panel. 
Several evolutionary tracks for different model masses are 
plotted. It is clearly visible that for a given metallicity, the 
masses of models which can enter the BL Her domain during 
the evolution are restricted. First, consider the lowest metal- 
licity models (Z = 0.0001, left-hand panel of Fig.rTTJ). Model 
with the largest mass displayed in the figure (M = 0.65 Mq) 
cross the BL Her domain at its lower part at L < 80 L©. 
Lower the model mass, higher the luminosity at the cross- 
ing with the BL Her domain. The models with M ^ 0.54 Mq 
enter the instability strip. For lower masses, crossing with 



4 Location of the instability strip is approximate. It was com- 
puted assuming M = 0.55 Mq, Z = 0.001 and convective param- 
eters of set PI. 



the instability strip occurs at larger luminosities, above the 
BL Her domain, or does not occur at all (evolution along al- 
most vertical tracks or even to the blue). As the metallicity 
is increased evolutionary tracks for models of a given mass 
(their pre- AGB part) shift towards the red and bend towards 
the lower luminosities. Consequently the less massive mod- 
els can enter the BL Her domain from the top. On the other 
hand, the tracks for most massive models leave the BL Her 
domain running below it. For Z = 0.001 (middle panel of 
Fig. [17| the tracks for models with 0.52 M Q sC M < 0.60 M 
cross the instability strip and for Z — 0.01 (right-hand panel 
of Fig. I17p the range of possible masses is even smaller, 
0.50 M < M s: 0.54 M . All tracks with M > 0.54 M Q 
run below the BL Her domain or even start beyond the red 
edge of the instability strip (0.60 M ). 

Qualitatively the same scenario is evident from the anal- 
ysis of the Dartmouth evolutionary tracks plotted in Fig. 1181 
The one-to-one comparison with the BaSTI tracks is diffi- 
cult, because the Dartmouth tracks were computed for spe- 
cific values of [Fe/H], which do not correspond to the Z 
values adopted in the BaSTI tracks. In Fig. [18] we have 
also plotted additional set of evolutionary tracks for mod- 
els adopting slightly enhanced abundance of the a-elements. 
Detailed list of the elements with their abundances can be 
found in table 2 of iDotter et alj <|200ct ). Their table 3 pro- 
vides the corresponding metallicity values, Z, which we also 
provide in Fig. 1181 Just as for the BaSTI tracks we ob- 
serve that lower the metallicity, larger the possible masses 
(and larger their range) of models crossing the BL Her do- 
main during the evolution. The main effect of increasing the 
abundance of the a-elements is shift of the tracks towards 
the lower luminosities and slightly lower effective temper- 
atures. Consequently, for the models with increased abun- 
dance of the a elements the range of possible masses that 
allow for crossing the BL Her domain (at given Z) shifts 
towards slightly smaller mass values. 

Pres ented results ar e also in qualitative agreement with 
the Y 2 (|Yi et al.l I1997T) evolutionar y models as well as 
with the Padova (|Girardi et alJl200Ch HB tracks (for M ^ 
0.55 Mq; tracks assuming lower M are not present in the 
Padova database). 

As discussed above, the evolution theory put constraints 
on the possible masses of the BL Her stars, which depend 
on the metallicity. Larger the metallicity, lower the possible 
masses of the BL Her stars. This has severe consequences 
for our best matching models for BLG184.7 133264 discussed 
in Section [3X2] A, B and C, all of which have Z = 0.01. 
The models are shown with different symbols in Figs 1171 
and !18l Comparison is straightforward for the BaSTI tracks, 
as these were computed also for Z = 0.01, and is done in 
the right-hand panel of Fig. 1171 We clearly see that mod- 
els B (M = 0.55 M ) and C (M = 0.60 M©) cannot be 
matched by the respective evolutionary tracks. Tracks for 
these masses run below the BL Her domain. To the contrary, 
model A is in excellent agreement with the evolutionary 
computations, the track for model of 0.50 Mq nearly crosses 
the square representing model A in the Figure. The one- 
to-one comparison with the Dartmouth tracks is somewhat 
difficult as metallicity of our best models, Z — 0.01, falls 
in between [Fe/H] = 0.0 and [Fe/H] = -0.5 for which the 
Dartmouth tracks were computed. Consequently, we have 
plotted the models in both the middle (Z ~ 0.0054) and 
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Figure 17. The BaSTI (Pictrinfcrn i et al ]|2006t) horizontal branch evolutionary tracks for three values of model's metallicity, Z = 0.0001 
in the left-hand panel, Z = 0.001 in the middle panel and Z = 0.01 in the right-hand panel. Evolutionary tracks start at the ZAHB, 
which is marked with the solid, horizontally running line. Each track is labelled with the corresponding model mass. Over-plotted are 
the edges of the instability strip, with the shaded area indicating the BL Her domain, with fundamental mode periods between 1 d (lower 
boundary of the shaded area) and 4d (upper boundary of the shaded area). 
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Figure 18. The Dartmouth llDotter et alj|2008tl horizontal branch evolutionary tracks for three values of model's iron abundance, 
[Fe/H] = —1.5 in the left-hand panel, [Fe/H] = —0.5 in the middle panel and [Fe/H] = 0.0 in the right-hand panel. The evolutionary 
tracks are labelled with the corresponding model mass. In each panel we also plotted additional set of evolutionary tracks for models 
assuming enhanced abundance of the a-elements, [ct/Fe] = +0.2 (dashed lines). Over-plotted are the edges of the instability strip, with 
the shaded area indicating the BL Her domain, with fundamental mode periods between 1 d (lower boundary of the shaded area) and 
4 d (upper boundary of the shaded area) . 



right-hand panel (Z ss 0.019) of Fig. [TU Metallicities are 
somewhat larger for evolutionary models assuming increased 
abundance of the a-elements (see labels in Fie. ll8[) . Conclu- 
sion is the same as for the BaSTI tracks. Only the least 
massive model A (M = 0.50 M©) can be matched with the 
evolutionary tracks. 

Luminosities of models B and C are much too high for 
their masses if confronted with the evolutionary tracks. Con- 
sequently, model A is our best model for BLG184.7 133264, 
selected based on both the pulsation computations and the 
evolution theory. 



5 DISCUSSION AND CONCLUSIONS 

In this paper we report the discovery of the first BL Her 
star showing the regular alternations of the light curve 



- the period doubling. The star was discovered nearly 
tw enty years after existence o f such objects was predicted 
bv lBuchlerfcMoskalikl da3), who found the period dou- 
bling behaviour in their radiative, BL Her-type models. 
The star, BLG184.7 133264, was found in the Galactic 
bulge OGLE-III survey data. It pulsates in the fundamen- 
tal mode, with period equal to 2.4 d. Its Fourier decompo- 
sition parameters firmly place it among BL Her variables. 
We also identified another strong period doubling candidate 
(BLG189.6 137529) in which however, the amplitude of al- 
ternations is very low and more accurate observations are 
necessary to confirm the detection. 

Discovery of period doubling in a BL Her star pro- 
vides a motivation for in-depth theoretical study of the phe- 
nomenon. In particular, the computation of new models, 
with up to date physics and including the convective energy 
transfer is of great value. In this paper we present the results 
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of initial model survey, which is focused on reproducing the 
period doubling behaviour in BLG184.7 133264. Thus, all 
computed models have Po fa 2.4 d and lie along an approx- 
imately straight line in the HR diagram. We have found a 
large domain, extending over a range of few tenths of solar 
luminosities, in which period doubling behaviour is possible, 
with rather large amplitude of alternations, of order of sev- 
eral hundreds of magnitude. The period ratios of the models 
showing the period doubling effect, indicate the crucial role 
of the 3:2 resonance between the fundamental mode and 
the first overtone in causing th e period doubling behaviou r, 
confirming the earlier result of iBuchler fc Moskalikl ([l992). 

The model light curves were transformed into the 
/-band and compared with the observed light curve of 
BLG184.7 133264. This was done with the intermediary of 
the Fourier decomposition parameters. We were able to 
match the observed pulsation amplitude and the amplitude 
of the period doubling only for relatively high metallicity, 
Z — 0.01. For three considered model masses, 0.50 Mq, 
0.55 Mq and 0.60 Mq, the model light curves reproduce the 
observed light curve quite well (Fig. I16|) . Low order Fourier 
parameters, R21 and $21 , are reproduced satisfactorily (par- 
ticularly for M = 0.60 M@). Also the secondary features of 
the light curve, in particular the pronounced bump at the 
minimum light are qualitatively reproduced. However, only 
the lowest mass model (M = 0.50 Mq) agree with the stel- 
lar evolution theory. For Z = 0.01 evolutionary tracks for 
models of M ^ 0.60 Mq run well beyond the red edge of the 
BL Her domain. Consequently, our best estimates for the 
parameters of BLG184.7 133264 are those of model A (Ta- 
bled}, specifically Z = 0.01, M = 0.50 M and L = 168 Lq. 

How robust are these estimates? Uncertainties lie in 
both the pulsation theory and in the evolution theory. Our 
conclusion about the high metallicity of BLG184.7 133264 
is based on the pulsation amplitudes of the computed mod- 
els. These depend on the adopted convective parameters, in 
particular on the mixing-length parameter, a, and the eddy 
viscosity parameter, a m . As we discussed in Section [3.3.21 
change of these parameters is not likely to affect our re- 
sults strongly. Nevertheless, more extensive model survey, 
including the variation of other parameters of the turbulent 
convection model, is needed and planned. In particular, the 
effect of including the turbulent pressure should be exam- 
ined. 

Derived metallicity can, in principle, be larger than Z = 
0.01, but only slightly. In test computations with Z = 0.02 
and assuming convective parameters of set PI, we find no 
traces of period doubling behaviour (independently of the 
model mass). 

At first glance, metallicity of order of Z = 0.01 may 
seem too high for the Population II object. The metal 
rich stars, however, are not exceptional in the Galactic 
bulge. Photometric metallicity estimates for several Type II 
Cepheids in the bulge range fro m —1.4 to +0.6, with mean 
value of [Fe/H] = -0.6 ± 0.17 l|Harris fc Wallersteinl [l98i 
IWallersteinll2002i ). Spectroscopic metallicity determinations, 
available mostly for bright K and M giants and Red Clump 
stars, indicate a large spread of metallicities in th e bulge, 
too, with [Fe/H] varying from —1.3 to +0.5 (see IZoccalil 
12011 for recent review). [Fe/H] determination for several 
main sequence stars in the bulge was also po ssible, dur- 
ing the microlensing events (jBensbv et alll201Ch . For these 



stars, metallicity range form —0.7 to +0.5. Metallicity can 
also be estimated for the RR Lyrae stars pulsating simul- 
taneously in the fundament al and the first overto ne modes 
(RRd stars). The OIII-CVS (|Soszvriski et al.ll201ll ) lists sev- 
eral bulge RRd stars with unusually short periods and low 
period ratios, down to Po ~ 0.35 d and Pi/ Pa ~ 0.726. The 
existence of these stars may be explained by high metal 
abundances, up to [Fe/H] w —0.36. As Z = 0.01 corresponds 
to [Fe/H] ~ —0.27, we see that our estimate of metallicity 
of BLG184.7 133264 is not unreasonable. 

For the luminosity, we can firmly indicate its lower 
limit. The period doubling domain in the HR diagram al- 
ways coincide with the crossing of the 3:2 resonant line with 
the Po = 2.4 d line. Run of the two lines shows very little 
sensitivity to the values of the convective parameters (Sec- 
tion l3.2[) . as well as to the mass and metallicity of the models 
(Fig. llO|l . The lowest luminosity at which the crossing occurs 
is w 150 L Q . 

Based on the recent evolutionary tracks, we have con- 
strained the mass of BLG184. 7 133264. With metallicity 
of Z = 0.01 only the models in a narrow mass range of 
~ 0.50 — 0.53 Mq evolve through the BL Her pulsation do- 
main. We get the best match with the pulsation models for 
M = 0.50 M . We note however, that at Z = 0.01 the evo- 
lutionary tracks that cross the BL Her domain vary rapidly 
with mass of the models. Although the current evolution- 
ary tracks of several groups lead to the same conclusions 
(Section [4} , the systematic errors in the micro-physics (e.g. 
opacities) or in the treatment of physical phenomena (e.g. 
convection, rotation, mass loss) may effect the results of 
evolutionary computations and somewhat change the best 
matching mass for BLG184.7 133264. 

Hydrodynamic model surveys (Bu chler fc Mo skalik 

this work) indicate that period doubling phenomenon 
should be a rather common property of the BL Her stars, at 
least o f those with periods close to 2.4 d. lBuchler fc Mo skalik 
( 1992) found the period doubling behaviour in a pe- 
riod range of 2.0 — 2.6 d. Predicted amplitudes of alter- 
nations are usually large, much larger than observed in 
BLG184. 7 133264. Thus, alternations should be easily de- 
tectable. Currently, around 25 BL Her variables are known; 
55 stars are listed in GCVS (ISamus et alj|2009h and close 
to 200 in OIII-CVS i|Soszvriski et al.l |200S| . |2010| . 2011, in 
prep.). OGLE observations are accurate enough to discover 
alternations with amplitudes as low as a few mmag. Hence, 
it is somewhat surprising that BLG184. 7 133264 is, at the 
moment, the only BL Her star with firmly detected period 
doubling behaviour. We note, that our strong period dou- 
bling candidate, BLG189.6 137529, also has a period in the 
expected range. On the other hand, only few of the known 
BL Her stars have periods close to 2.4 d. Our study shows 
that higher metallicity (e.g. close to solar, which may not 
be exceptional for BL Her stars) may inhibit the period 
doubling effect. Definitely, more detailed non-linear study 
is needed. 

In the forthcoming paper we plan to study the full grid 
of models, covering the whole BL Her domain in the HR dia- 
gram, not only the 2.4 d models. Such survey will reveal the 
structure of the period doubling domain in the HR diagram, 
in particular the periods at which the phenomenon is possi- 
ble and factors which control the existence and strength of 
the period doubling effect. 
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